Parameterization dependence of T matrix poles and eigenphases 
from a fit to nN elastic scattering data 
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We compare fits to -kN elastic scattering data, based on a Chew-Mandelstam Tf-matrix formalism. 
Resonances, characterized by T-matrix poles, are compared in fits generated with and without 
explicit Chew-Mandelstam Ji'-matrix poles. Diagonalization of the S matrix yields the eigenphase 
representation. While the eigenphases can vary significantly for the diflferent parameterizations, the 
locations of most T-matrix poles are relatively stable. 
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I. INTRODUCTION 

The excited states of the nucleon [l| have been stud- 
ied in a wide array of reactions initiated mainly by pion 
and photon beams. Other approaches have involved an 
examination of the invariant mass distribution ofprod- 
ucts from, for example, nucleon-nucleon reactions J3] and 
J/* decays [3]. Most states listed by the PDG [if were 
identified from fits to -kN elastic scattering and reaction 
data. Photo-decay amplitudes were determined mostly 
through analyses of single-pion photoproduction data. 

Recent measurements of cross section and polarization 
quantities, related to the photo- and electroproduction 
of states other than ttN , have been analyzed separately 
and in multi-channel approaches. These studies have pro- 
vided stronger evidence for states seen only weakly in -kN 
elastic scattering, and have suggested new states, cou- 
pling more strongly to other channels 

Among the most extensive ttN scattering analyses [H- 
0] , the parametrization of Ref . Q based on the SAID in- 
teractive fitting and database codesQ (the SAID-GW 
fit), utilizing the most recent data, has found the fewest 
number of N and A resonances. In the fit of Ref. [9,], a 
search for weaker structures was carried out. There, the 
existing solution was modified using a simple product S- 
matrix approach, to include the effect of an added Brcit- 
Wigner resonance in each partial wave. Chi-squared was 
mapped for various combination of masses, widths and 
branching fraction. Two marginally significant candi- 
dates were found in the and i^is partial waves, with 
pole positions: 1689 - i96 MeV (for Su) and 1793 - i94 
(for Fis). Of these, the i^i5 has been reported in subse- 
quent fits, while the has not. 



Here we have considered another approach. As de- 
tailed further below, existing GW-SAID fits to ttA^ elas- 
tic scattering data have utilized a fit form based on the 
Chew-Mandelstam (CM) if-matrix. This approach is ca- 
pable of generating T-matrix poles without the assump- 
tion of explicit fC-matrix poles. Previous fits have only 
included an explicit CM i^-matrix pole for the A(1232). 
In the present study, an alternative parametrization with 
one explicit CM X-matrix pole in each partial wave 
is used to generate a fit independent of the usual CM 
parametrization form. 

The motivation for this new fit is twofold. By changing 
the parameterization, we are able to gauge the stability 
of the amplitudes and resonance positions. We are also 
able to see if the addition of new explicit CM if-matrix 
poles translates into additional resonance signals. 

Below, in Sec. |lTl we briefly review the CM _fC-matrix 
formalism used in this and previous fits. The eigenphase 
representation, and some numerical details, are reviewed 
in Sec, mil Results for the partial wave fits and resonance 
spectrum are compared in Sec. IIVI Finally, in Sec. |Vl we 
consider the implication of this and future work. 



II. CHEW-MANDELSTAM FORMALISM 

The Chew-Mandelstam approach for the parametriza- 
tion of multichannel hadronic ttN elastic scattering and 
reactions to other hadronic channels has been described 
in detail in Refs. The x^-fits to data have been 

additionally constrained using the forward dispersion 
relations and fixed-t dispersion relations for the invariant 
B amplitudes. 

The standard CM parametrization can be expressed in 
terms of the on-shell Heitler partial wave K matrix, K 
as 



K- 



\E) = K ^{E)-KcC{E), 



(1) 
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where E is the (complex) scattering energy, K is the 
Chew-Mandelstam K matrix and C is a diagonal matrix, 
whose matrix elements are termed the Chew-Mandelstam 
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functions [l3|. The Heitler K matrix is related to the 
partial wave transition amplitude matrix, T as 



T-\E)^K-\E)-ip{E). 



(2) 



Here, p{E) = 5{E—Hq), where is the (model indepen- 
dent) relativistic free-particle Hamiltonian with physical 
(stable) particle masses. It determines the CM functions, 
C{E) via the relation lmC{E) = p(S)0 

The standard form used in the GW fits is defined by 
the choice for the CM if-matrix elements 



if(ii;) = ^c„^"(i?), 



(3) 



where Cn are a set of constants and z is a linear function 
of the scattering energy, E. The integer, n is typically 
between 2 and 5, and depends on the matrix element in 
question. 

Note that K defines an entire function of the complex 
parameter E for finite values. This form is used for all but 
the P33 partial wave, which includes an explicit pole in K. 
For partial waves other than the P33 , we see that the CM 
K matrix, K, is without poles (or other singularities). 
The Heitler K matrix, K, 



K 



1 



1 - A'[ReC] 



-K 



(4) 



has a pole whenever det[l - YleC{E)K{E)] = 0. The 
matrices K and K are free of branch point singulari- 
ties 

The alternate form of the CM K matrix is similar to 
the form used in the P33 partial wave of the standard K 
parametrization, described above. This form is given by 



K, 



1^1 

E — En 



(5) 



Here, ^i{E) is a polynomial without a zero at the pole 
position, _Ep, and the index labels the channel {t^N , ttA, 
pN , and 'qN), fi{E) is an entire function of the complex 
energy, E. 



III. EIGENPHASE REPRESENTATION 

The fit produces a unitary S-matrix of amplitudes for 
all contributing channels. While those channels not fitted 
to data are unlikely to give a quantitative representation 
of the reaction (for example ttN — j> ttA), they can be 
used to construct a set of eigenphases, which provide an 
interesting characterization of resonance behavior. 



^ The inclusion of quasi-two-body channels, such as vrA are con- 
strained to have branch points at the stable three-body thresh- 
olds. 



The unitarity of the 5* matrix implies that its eigen- 
values are phase factors. The matrix, U , of eigenvectors 
diagonalizes the S matrix as 



U'^SU = A, 



where 



A = 



/Ai 
A2 












(6) 



(7) 



Exploiting ~ 1 we write 



A,, 



(8) 



with (pi real. 

Our objective is the numerical evaluation of the eigen- 
phases given the T-matrix elements from various fits. 
This is straightforward at a given energy, using a stan- 
dard routine to diagonalize the unitary S matrix. The 
only complicating issue is correlating a given eigen- 
phase, 4'i{E) with the appropriate eigenchannel when two 
(or more) eigenphases converge as the energy changes. 
In other words, once an eigenchannel i is determined, 
we must track it for all energies. The no-crossing 
theorem [3| is readily generalized to unitary matrices and 
shows that, in a given partial wave, the eigenphases may 
not be equal for any energy. This property is exhibited 
in the eigenphase plots discussed below. 

Given the T matrix at some energy, T{E), we can form 
the S{E) matrix. We diagonalize this matrix using a 
standard routine to obtain the eigenvalues {Xi{E)}2^^, 
where n is the number of channels. 

If the eigenvalues are nearly degenerate at some energy, 
it is difficult to distinguish which eigenvalue corresponds 
to a given eigenchannel, say i, since diagonalization of 
S doesn't preserve the eigenchannel ordering. The set of 
eigenvectors, however, must be orthogonal at any energy; 
and, for continuous partial wave amplitudes, the change 
of the eigenvector for a given eigenchannel is small for 
nearby energies. 

The eigenchannels are maintained using the following 
method. The S matrix is diagonalized at the initial en- 
ergy, say El = 1150 MeV. We obtain n eigenphases 
(where n is the number of channels included for the given 
partial wave), Ai(i?i), . . . , A„(i?i) and their correspond- 
ing eigenvectors vi{Ei), . . . ,Vn{Ei). We wish to correlate 
the eigenvalues and eigenvectors with a given eigenchan- 
nel throughout the evaluation of the eigenvalues at higher 
energies, E > Ei. 

Increasing the energy a small amount (10 — 15 MeV) to 
E2, we again diagonalize the S matrix and evaluate the 
Xi{E2), Xn{E2) and eigenvectors wi(£^2), • ■ . ,w„(£^2)- 

In order to track the eigenchannel, we evaluate the 
matrix of overlaps: 



0,j{EuE2) = v,{Ei)^Vj{E2). 



(9) 
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FIG. 1. (Color online) Selected partial-wave amplitudes (L2/,2j)- Solid (dashed) curves give the real (imaginary) parts of 
amplitudes corresponding to the WI08 Q solution. Dash-dotted (dotted) curves give the real (imaginary) parts of amplitudes 
corresponding to the XP08 solution, (a) Sn, (b) S31, (c) Pn, (d) -D13, (e) F15, and (f) F37. All amplitudes are dimensionless. 
Vertical arrows indicate resonance Wr values and horizontal bars show full F and partial widths for F.^jv. 



As E2 El, we have 

lim O,, {El, E2)^6,j, (10) 

which is just the statement that the eigenvectors are or- 
thonormal. For E2 — Ei ~ 10 MeV, we identify the 
eigvenvalues according to the largest overlap in the set 

m,iEi,E2)\}]^i. (11) 

Suppose, for example, that we have three channels and 
at the energy Ei, we write the eigenvalues in the order: 

Ai,A2,A3. (12) 



And at energy E2 for i — I, we find that 

\Oi3{Ei,E2)\ > \OniEi,E2)\ > |0i2(Si, -^2)!, (13) 

then for energy E2 we order the eigenvalue A3 first; the 
ordering for the other eigenvalues is determined similarly. 

IV. RESULTS 

The fits with (XP08) and without (WI08) exphcit CM 
/C-matrix poles, in waves other than (^2/^27) -P33J are 
compared in Fig. [TJ Differences in the partial waves are 
slight, and the fit quality is comparable over the reso- 
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FIG. 2. (Color online) Eigenphases. (a) Sn, (b) S31, (c) Pn, (d) D13, (e) F15, and (f) Fg' 



nance region, each fit using a similar number of param- 
eters. This feature of the ttN elastic scattering analysis 
seems quite stable. 

In Fig. [21 we have calculated the eigenphases corre- 
sponding to the full S-matrix. Only the ttN ttN and 
ttN TjN channels have been constrained by data. The 
behavior of these phases does, however, provide an inter- 
esting perspective on the emergence of resonance struc- 
tures in the fits. 

In the Sii partial wave, both fits have two eigenphases 
crossing 90 degrees, at 1535 and 1670 MeV for WI08, 



and at 1530 and 1700 MeV for XP08. If one computes 
the usual Heitler K matrix, as was done in Ref. [Tg} , 
ii'-matrix poles are found at these energies (since the 
unitary transformation, U \Ei(\.^] diagonalizes K simul- 
taneously with S and Ku — tant/)^). In the partial 
wave, a 2- and 3-channel fit are compared, yielding identi- 
cal crossing energies, again corresponding to a K-vaaXrvx. 
pole (at 1655 MeV). Note that in the WI08 plot, two 
eigenphase nearly touch, but do not cross. 

In the Pn plot, only one of the solutions has a 90 de- 
gree crossing leading to a if-matrix pole. Note, however. 
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TABLE I. Pole positions in complex energy plane of the T matrix for the nN — >■ -kN reaction. The functional forms (see text) 
employed in the SAID fits are compared for selected partial waves. Each T pole position is expressed in terms of its real and 
imaginary parts [Mr, —Tfi/2) in MeV. The second sheet pole is labeled by a f. 
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(1647,42) 
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that the energy dependence of the eigenphase crossing, 
and nearly crossing 90 degrees, is very similar. This fea- 
ture determines another measure of resonance behavior, 
to be discussed below. 

The Di3 eigenphases are quite different in the two fits. 
In the WI08 fit, there are no 90 degree crossings, while in 
XP08, we see two crossings. This hints at a different res- 
onance structure, though the nN T matrices are nearly 
identical. 

In the Fi5 and F^y eigenphase plots, the XP08 solution 
has a single crossing, whereas the WI08 solution does not. 
Here also, a comparison of the eigenphases which cross, 
or come close to crossing, 90 degrees have a similar energy 
dependence. 

As has been noted previously resonances may be 
associated with a single eigenphase crossing 90 degrees, 
and this will result in a i4'-matrix pole. However, a more 
robust measure (if a set of amplitudes is available) is 
given by the time-delay matrix [isj . which is propor- 
tional to the sum of energy derivatives of all eigenphases. 
Other factors, such as threshold openings can also pro- 
duce rapid energy dependence. Certainly the correct 
method of resonance identification requires the location 
of poles in the complex energy plane on unphysical sheets 
close to the physical region, which we demonstate be- 
low. Our employment of the eigenphase approach illus- 
trates the fact that resonance structure may vary with 



the appearance of resonances in different parametriza- 
tions without significantly altering the shape of the t:N 
elastic amplitude. It is usually the case, however, that 
such resonances are deep in the complex plane having 
large widths. 

In Fig.[3l for illustration, we plot the sum of eigenphase 
energy derivatives for the Pu and -D13. The peaks for 
Pii are nearly identical and occur at about 1350 MeV, 
which (we will see) corresponds with the real part of the 
pole position. For the D13, peaks corresponding the the 
4-star state, near 1500 MeV, are closely aligned. The 
second peak has almost no evidence in the ttN elastic 
amplitude. However, a large contribution to the (unfitted 
and therefore unconstrained) ttA or pN channels would 
result in the second peak. 

In Table U we compare the pole positions associated 
with resonance behavior in the plotted amplitudes. The 
third 5*11 pole in XP08 closely resembles the structure 
found in Ref. at (1689,96) MeV, by scanning all par- 
tial waves with an added Breit-Wigner contribution. The 
very broad (1646, 290) MeV Pu state is similarly close 
to one found in the SM90 fit at (1636,272) MeV. 
Two extra poles were found in the D13 partial wave 
for the XP08 solution compared to WI08. We do not 
intend to report the (1716, 370) MeV pole as a reso- 
nance but merely mention it here in connection with 
the present sensitivity study. Interestingly, the pole at 
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(1740, 66) MeV has its effect masked by a zero interven- 
ing between the pole and real energy axis and therefore 
makes little impact in the physical region. 

V. CONCLUSION 

We have reported a study of the parameterization de- 
pendence of our -kN elastic amplitudes and resonance 
spectrum, using a very different form for the CM K ma- 
trix, with explicit poles in each partial wave. The partial- 
wave amplitudes were found to be very stable with the 
change. 

The eigenphase representation was introduced, mainly 
as a novel approach to resonance identification, and be- 
cause it provides a more concrete example of properties 
discussed in older works. This discussion also provides a 
continuation of the study started in Ref. [T6|. 

The more formally correct extraction of pole positions 



has revealed structures mainly found in earlier fits to the 
■kN elastic scattering data. As the partial wave ampli- 
tudes have not changed significantly, the effects of new 
resonances must be minimized through large widths, in- 
tervening zeros, or small coupling to the TriV channel. 

In the SM90 fit, a study of the resonance spectrum was 
tried where, in addition to experimental data, the ampli- 
tudes from the KH 5] and CMB [6] analyses were added 
as soft constraints. A possible extension to the present 
work would be a re-examination of the resonance spec- 
trum from a fit, with explicit CM iC-matrix poles, con- 
strained to more closely follow either the KH and CMB 
analyses, or a multi-channel analysis. 
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